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We discuss the region of transition between galactic and extragalactic cosmic rays. The exact shapes and 
compositions of these two components contains information about important parameters of powerful astrophysical 
sources and the conditions in extragalactic space. Several types of experimental data, including the exact shape 
of the ultrahigh energy cosmic rays, their chemical composition and their anisotropy, and the fluxes of cosmogenic 
neutrinos have to be included in the solution of this problem. 



1. INTRODUCTION 

We believe that at some very high energy of 
order lO-'^^-lO^^ eV all observed cosmic rays come 
from extragalactic sources because they cannot 
be contained in the Galaxy long enough to be 
accelerated to such high energies [1]. The cur- 
rent knowledge of the cosmic ray spectrum and 
composition do not allow us to determine where 
exactly the transition from local (galactic) to ex- 
tragalactic origin happens. 

We know that extragalactic cosmic rays lose 
energy in propagation from their sources to us if 
their sources are isotropically and homogeneously 
distributed in the Universe. The main energy loss 
process is the photoproduction interaction in the 
microwave background (MBR) that causes the 
GZK effect, the steepening of the cosmic ray spec- 
trum above 6x10^^ eV [5]. Several fits of the ex- 
isting experimental data have been published in 
recent years that derived different values of the 
most important astrophysical and cosmological 
parameters: the acceleration (injection in terms 
of cosmic ray propagation) energy spectrum of 
these particles and the maximum acceleration en- 
ergy, their chemical composition and the cosmo- 
logical evolution of their astrophysical sources. In 
the assumption that most of ultrahigh energy cos- 
mic rays (UHECR) are protons injection spectra 
as different as E"^ '' [3] and E"^ " [4] and cos- 
mological evolution of the type (1 -I- z)"* with m 
values from to 3 have been obtained. 



Other attempts |5|6|7j have assumed that ex- 
tragalactic cosmic rays have at their sources the 
same mixed chemical composition as the low en- 
ergy galactic cosmic rays. In such a case the 
main energy loss process is the disintegration of 
the heavy nuclei mostly in interactions in MBR. 
Hadronic interactions become important only af- 
ter the energy per nucleon exceeds the photopro- 
duction threshold. The fits of the observed cos- 
mic rays spectrum under this assumption gives 
an intermediate £;-2-2-2.4 jj^jg^tion spectrum. 

In all these attempts the fits of the observa- 
tions show the end of the galactic cosmic rays 
spectrum which is obtained by subtraction of the 
propagated extragalactic spectrum from the ex- 
perimentally observed one. This process gives 
some limits of the astrophysical parameters [8] 
when the subtraction gives unphysical negative 
values. 

In addition to the cosmic ray spectrum there 
are several types of relevant data: UHECR chem- 
ical composition, the anisotropy in this energy 
range, and the production of cosmogenic [9J neu- 
trinos by extragalactic cosmic rays. We will 
briefly discuss these parameters. 

2. UHECR ENERGY SPECTRA AND 
COMPOSITION 

Figure [T] compares two different fits of the ex- 
tragalactic cosmic ray spectra in the assump- 
tion that they are purely protons and that 
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their differential acceleration spectrum is a power 
law E~'^~^^. The experimental data are from 
AGASA [TU] and HiRes [TT] and are normalized 
to each other at 10^^ eV. Since we are now only 
interested in the shape of the spectrum the exact 
differential flux at the normalization points is not 
important. It is obvious that the two experimen- 
tal measurements agree well with each other on 
the shape of the energy spectrum with exception 
of the AGASA events above 10^° eV. The most re- 
cent analysis of the AGASA data, presented this 
summer [12] decreases the energy assignment of 
the AGASA data by 10-15% and makes the spec- 
trum closer to that of HiRes. 

Fit a |3] derives an injection spectrum with 7 = 
1.7. The dip at about 10^^ eV is due to the transi- 
tion of proton energy loss to Bethe-Heitler e+e~ 
pair production to purely adiabatic loss as pre- 
dicted by Berezinsky&Grigorieva [13^ . The model 
does not need any contribution from galactic cos- 
mic rays to describe the observed cosmic ray spec- 
trum down to 10^^ eV. There is also no need for 
cosmological evolution of the extragalactic cosmic 
ray sources, although some source evolution can 
be accommodated with a slight change of 7. The 
model predicts purely proton composition of the 
extragalactic cosmic rays and does not work as 
well as shown in Fig. [1] if more than about 10% 
of the cosmic rays at the source are nuclei heav- 
ier than H. Fit b j4] does need contribution from 
galactic sources with i<;~3-50 spectrum that ex- 
tends well above 10^^ eV. The extragalactic con- 
tribution is shown for two different cosmological 
evolutions with m = 3 (as used in Ref. [5) and 
4 (upper edge of the shaded area). The influ- 
ence of the cosmological evolution on the cosmic 
ray propagation is modest because no protons 
injected at redshifts z higher than 0.4 arrive at 
Earth with energy above 10^^ eV independently 
of their initial energy. 

Obviously these two models predict very dif- 
ferently the end of the galactic cosmic ray spec- 
trum. In model a the galactic cosmic ray sources 
do not need to accelerate particles above 10^* eV. 
In model b they should be able to reach energies 
higher by one and a half orders of magnitude. 
This would affect very strongly the expected cos- 
mic ray composition. 
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Figure 1. Comparison of two fits of the UHECR 
spectrum. The solid line shows model a from 
Ref. [3J and the dashed lines show model b from 
Ref. [4J. Short dashes show the galactic compo- 
nent needed to fit the observed spectrum in b. 

The assumption that extragalactic cosmic rays 
have a mixed composition at acceleration gives 
somewhat intermediate results for the injection 
spectrum of UHECR [5|6|7j . The spectrum that 
fits the observation best has 7 values between 1.2 
and 1.4. The chemical composition of cosmic rays 
at Earth is also different and obviously depends 
on the source composition. 

It used to be that we considered three sepa- 
rate regions in the cosmic ray spectrum - energy 
below 10^^ eV where the acceleration is at su- 
pernova remnants, 10^^-10^^ eV where acceler- 
ation is at unknown galactic sources and extra- 
galactic component. Since stochastic shock ac- 
celeration is rigidity dependent (as is the escape 
from the Galaxy) this picture predicts a compli- 
cated composition energy dependence. When a 
source cannot accelerate protons any more, there 
are only heavier nuclei, initially He, then CNO, 
then Si and finally Fe nuclei that are acceler- 
ated. One may imagine the composition becom- 
ing heavier between 10^^ and 10^^ eV, then lighter 
at higher energy when the unknown sources take 
over, then heavier again when they are exhausted, 
then lighter again when extragalactic cosmic rays 
start dominating. Recent developments in the ac- 
celeration theory [14j show that the maximum ac- 
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celeration energy may be much higher and 
galactic sources can cover the whole energy range 
to 10^^-10^^ eV, i.e. we expect one single transi- 
tion region. The Kascade experiment has already 
observed [16] the beginning of the exhaustion of 
the galactic cosmic ray sources and the respective 
change of the composition. 

Figure [2] shows the predicted composition as a 
function of the total energy per nucleus in the 
transition region in the classical cosmic ray units 
of < In A >. This is very appropriate when the 
detection is through air showers with logarith- 
mic sensitivity to both energy and composition. 
Model a presents the easiest case to explain. All 
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Figure 2. Cosmic ray chemical composition as 
a function of the total energy per nucleus in the 
three models. The mixed composition model data 
are taken Ref. [6]. The error bars assumed are 
logarithmic and probably lower than the true ex- 
perimental ones. The very heavy dashed line rep- 
resents symbolically the results of the HiRes pro- 
totype/MIA coincidences. 

composition changes happen below 10^^ eV. We 
do not plot the composition below that energy 
as it is not very well defined. The results for the 
other two models are presented using our own un- 
derstanding of them. For model b we assume that 
at 10^^'^ eV all cosmic rays are iron and assign 
< InA > value of 4. At higher energy we use 
the fraction of galactic cosmic rays to iron nuclei 
and thus calculate the corresponding < InA > 
value. For the mixed composition we use Fig. 3 



of Ref. [6^ where the cosmic ray composition at 
Earth is plotted for 7=1.3 and E^^^ is given as 
1020-5 eV. 

Surprisingly the difference in the energy depen- 
dence of < InA > of model b and of the mixed 
composition model is not that big, although in 
the model b case the composition only consists 
of a combination of Fe and H, and in the mixed 
composition case we have five groups of nuclei. 

Model a gives a very different picture, at least 
above 10^^ eV, where the composition is purely 
Hydrogen. It is definitely distinguishable from 
the other two models. The prediction that the 
composition does not change above 10^* eV and 
is very light is supported by the HiRes X^ax 
measurement [17J, indicated with heavy dashes 
in Fig. [21 What is the exact meaning of the de- 
tected light composition is not known because of 
the differences between the hadronic interaction 
models used for data analysis. 

On the other hand, other experiments support 
a much milder energy dependence of the cosmic 
ray chemical composition, that is more in line 
with the prediction of of models b and that of 
mixed extragalactic cosmic ray composition. 

3. ANISOTROPY 

Low energy cosmic rays diffuse in the magnetic 
fields of the Galaxy and lose memory of the lo- 
cation of their sources. The anisotropics are very 
small, well below 1% and are difficult to measure. 
The measurement of a small anisotropy with air 
shower experiments requires an exact knowledge 
of the experimental acceptance and of the lifetime 
of the shower array. 

At higher energy things are supposed to change 
as the cosmic ray rigidity increases and their dif- 
fusion in the Galaxy becomes faster. The first 
question we still cannot answer is of the rigidity 
(energy) dependence of the cosmic ray diffusion 
coefficient. An energy dependence of £'~0-5-0-6 
is derived from the ratio of secondary to primary 
cosmic rays. If this dependence is extended by 
seven orders of magnitude then galactic UHE pro- 
tons should show very strong anisotropy that is 
not observed. The two spots with increased cos- 
mic ray flux identified by the AGASA group [TS] 
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(one in the general direction of the Galactic cen- 
ter and one in the direction of Orion) are not 
confirmed by other experiments. Theoretically 
we expect Kolmogorov turbulence in the Galaxy, 
that should under some circumstances give E^^^ 
energy dependence. 

One way to explain the low anisotropy in the 
beginning of the transition region is to assume 
that all galactic cosmic rays are heavy nuclei. If 
they were indeed heavy nuclei the average parti- 
cle rigidity would significantly decrease and would 
maintain the high isotropy of the galactic cos- 
mic rays. In the transition region, however, the 
cosmic ray chemical composition becomes lighter 
and correspondingly we expect to see some degree 
of anisotropy. 

I am convinced that in the case of high exper- 
imental statistics (as expected from the Pierre 
Auger Observatory [19]) the anisotropy in the 
transition region between galactic and extragalac- 
tic cosmic rays deserves a careful study. A part 
of it is theoretical. We should propagate particles 
of different rigidity in the Galaxy and attempt to 
understand their general behavior. Analytic so- 
lutions of the diffusion equation are not any more 
suitable for this problem. Particle trajectory has 
to be numerically solved, most likely in a Monte 
Carlo fashion, and in detailed enough magnetic 
field models. Such models should be tested to 
match analytic calculations when applied to the 
appropriate simplified models. 

4. COSMOGENIC NEUTRINOS 

Cosmogenic neutrinos are produced in the same 
photoproduction interactions of the UHECR pro- 
tons that create the GZK effect. They were 
first proposed by Berezinsky&Zatsepin in 1969 [9 
and were the subject of many calculations after- 
ward. Cosmogenic neutrinos are often considered 
a 'guaranteed source' of ultrahigh energy neutri- 
nos. They are indeed guaranteed, since we know 
UHECR exist, but their flux is uncertain. 

An essential quality of neutrinos is that they 
have a low interaction cross section. This makes 
neutrino detection a difficult problem that re- 
quires huge detectors of at least km^ scale. Such 
detectors are now in the stage of planning [5D] and 



construction [21|22j . The question of the cosmo- 
genic neutrino fiux and of its relation to other 
data from ultra-high astrophysics experiment is 
very timely. 

The low interaction cross section is not only 
a deficiency. While protons emitted at redshifts 
z exceeding 0.4 do not reach us with energy 
above 10^^ eV the cosmogenic neutrino produc- 
tion peaks at redshifts exceeding 2 for (1 + z)^ 
cosmological evolution of the UHECR sources. 
This is the main link between the extragalactic 
UHECR and cosmogenic neutrinos. The detec- 
tion of cosmogenic neutrinos may help the degen- 
eracy in modeling of the extragalactic cosmic rays 
spectra shown in Fig. [T] 

The point is that models that use flat cosmic 
ray injection spectrum and require strong cosmo- 
logical evolution of the UHECR sources, such as 
model h would produce significantly more cos- 
mogenic neutrinos than steep injection spectrum 
models with no cosmological evolution. Such 
a comparison between the two models is shown 
in Figure [3l The figure shows cosmogenic neu- 
trino fiuxes generated by proton interactions in 
MBR [53] and in the infrared/optical background 
(IRB) [23] with the injection spectrum and cos- 
mological evolution of the two models as indi- 
cated in the two panels. The injection spectra 
of the two models are normalized to the UHECR 
differential flux at 10^^ eV. 

Before discussing the magnitude of the two 
fiuxes I want to attract your attention to the flux 
of electron neutrinos and antineutrinos that ex- 
hibits two maxima separated by about two orders 
of magnitude of energy. The higher energy one is 
due to the i:'^ i^fj, + i^fi + i^e decay final state and 
consists mostly of electron neutrinos. It peaks at 
the same energy where fi^ and do. The lower 
energy maximum is due to electron antineutri- 
nos from neutron decay. Comparison between the 
neutron interaction length and their decay length 
show that all neutrons of energy below 10^°'^ eV 
decay rather than interact. 

The number density of IRB is of order 1 and 
varies by about a factor of 2 in different mod- 
els, but its energy extends to energies exceeding 
1 eV. What that means is that extragalactic par- 
ticles of much lower energy would interact in IRB 
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Figure 3. Fluxes of cosmogenic neutrinos gener- 
ated in proton interactions in the isotropic photon 
background [2i[ in the two models of the UHECR 
spectrum. The upper panel compares the fluxes 
of + Pg. The lower one compares the fluxes of 



and will generate neutrinos. Reference [24] shows 
the neutrino yields generated by protons of en- 
ergy as low as 10^^ eV. This yield is small but 
has to be weighted by the much higher number 
of particles at that energy - almost 1,000 higher 
that that of 3x10^^ eV for a flat E^^ injection 
spectrum. A more recent calculation that also 
includes UV photons [27] employs interactions of 
10^^ eV nuclcons that have to be weighted still 
higher. The cosmological evolution of IRB is not 
as strong as that of MBR, as the radiation is gen- 
erated through emission directly by stars and af- 
ter absorption by dust, but still significant up to 
redshift of 3. The IRB model of Ref. [21] is used 
for the shown calculation. 

Model b generates much higher cosmogenic 



neutrino fluxes than model a because of two rea- 
sons that contribute roughly the same increase of 
the cosmological neutrinos. Firstly, it uses much 
flatter injection spectrum E~'^''^ which means 
equal amount of energy per decade. It thus 
contains many more particles above the photo- 
production interaction threshold in the MBR is 
about 3x10^^ eV and, of course, decreases as 
il + z)-\ 

The other reason is that model b employs a 
strong cosmological evolution of the cosmic ray 
sources. This increases by ^3 the number of par- 
ticles injected at redshift of 2, but in addition 
increases by a factor of three the number of par- 
ticles above the interaction threshold. This way 
the total neutron flux is increased by the cosmo- 
logical evolution of the sources by a factor of five. 

The difference in the peak values of the cos- 
mological neutrinos generated by the two models 
is more than one order of magnitude. In practi- 
cal terms, after the energy dependence of the i^N 
cross section is taken into account, this means 
that model b generates fluxes that are in princi- 
ple detectable by IceCube at the rate of roughly 
less than one event per year, while model a gener- 
ates undetectable fluxes of cosmogenic neutrinos 
in km'^ detectors. 

Ice and water neutrino detectors are generally 
not very suitable for cosmogenic neutrino detec- 
tion. Much better strategies for these UHE neu- 
trinos are the radio and acoustic detectors that 
have very high detection threshold, but also will 
have higher effective volume. The other option 
are giant air shower arrays such as Auger, that 
can reach effective volume of 30 km"^ and, with 
sufficiently low threshold (10^^ eV) could detect 
several events per year. 

Cosmogenic neutrinos are also generated in the 
mixed composition scenario [25|26j . Since the 
major energy loss process is the disintegration of 
heavy nuclei, the main flux component of 9^ from 
neutron decay. The D^. flux exceeds by a factor 
of 5 the sum of the neutrino fluxes of all other 
flavors in Ref. [26]. The absolute magnitude de- 
pends again mainly on the injection spectrum and 
the cosmological evolution of the sources. 
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5. DISCUSSION 

The transition region between galactic and ex- 
tragalactic cosmic rays is currently not very well 
defined and studied. It has a significant astro- 
physical importance, because the energy content 
of the extragalactic cosmic rays, as well as their 
injection spectrum and composition, will define 
much better the type of their sources. In addi- 
tion, a derivation of the cosmological evolution 
of the sources will not only restrict the number of 
available scenarios, but also provide an additional 
measure of the cosmological evolution of powerful 
astrophysical objects. 

The first experimental result on the shape 
of this transition presented by the Fly's Eye 
group [29] used the simultaneous change of the 
cosmic ray composition and energy spectrum 
shape. We should now employ all available types 
of measurements to solve this problem. These 
measurements include the cosmic ray composition 
in the transition region and above it, the cosmic 
ray anisotropics, and the possible detection of cos- 
mogenic neutrinos. 

The possible detection of cosmogenic neutrinos 
would be a powerful test if we succeed in collect- 
ing a reasonable statistics of such events. It is un- 
likely this will happen very soon, but the expand- 
ing efforts for designing and building detectors for 
ultrahigh energy neutrinos are very encouraging. 
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